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Py ou Ce 


The TACTICS computer presram developed by the RAND 
Corporation has the capability of simulating three-body, three- 
dimensional tactical situations. This study investigates the problem 
of using the TACTICS program when one of the bodies is to have the 
characteristics of a helicopter. The differences between the forces 
acting ona helicopter and the forces assumed by the TACTICS 
program are discussed, and a method to Seats the helicopter ° 
is presented. The basic idea of this method is to ''fly"' a hypo- 
thetical fixed-wing aircraft with parameters which enable it to 
perform like a helicopter. The programing required to simulate 
several basic helicopter maneuvers is discussed, and two examples 


of a helicopter ina tactical environment are presented. 
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I, INTRODUCTION 


ie BACKGROUND 

A three-body, three-dimensional simulation program entitled 
TACTICS was developed at the RAND Corporation primarily as an 
aid in the study of the close-in phase of two maneuvering fighter air- 
craft in a combat situation. The program was constructed ina 
modular form consisting of many subroutines, allowing the user to 
make modifications to the program for other areas of study. Reference 
] indicates possible applications of the program to the simulation of 
air-to-air, surface-to-air, and air-to-surface missiles, as wellas 
ICBM reentry vehicle, rocket booster, and satellite investigation. 

In recent years, the employment of the helicopter in conventional 
Petare bas increased significantly. The missions of the helicopter 
now include not only the familiar search and rescue roles, butas 
mentioned in Ref. 2, profiles that cover the air-to-ground environ- 
ment, including suppressive fire at area, light, and hard point targets. 
Helicopters have become fast and highly maneuverable, some having 
the capability of self-protection in hostile air and ground battle situa- 
tions. While confrontations between helicopters and opposing force 
air-to-air and major surface-to-air defenses have been rare, in- 
cidents of this type may be expected to continue, and possibly 
increase, in future limited conflicts. 

An analytical study of the helicopter ina tactical situation 
would attempt to provide answers to such diverse problems as 


optimum tactics to be employed by both the helicopter and opposing 








weapon system, missile and sensor design parameters, and trade- 
offs between helicopter maneuverability and defensive armament. 
The TACTICS program has the capability of providing an output with 
vehicle position, velocity, acceleration, altitude, geometry, and 
aerodynamics listed by time history. However, in order to be able 
to provide accurate output, the program must be initialized with data 
that correctly describes the characteristics of the vehicle under con- 
sideration. These characteristic parameters have been developed 
quite thoroughly for certain high-speed vehicles, but the TACTICS 
manual, Ref. 1, does not indicate application of the program to 


studies of a vehicle with the unique capabilities of a helicopter. 


ae OBJECTIVE 

It is the objective of this study to devise a method for applying 
the TACTICS program to the problem of simulating a helicopter ina 
tactical environment, such as in an encounter witha high perfor- 


mance aircraft or when under attack froma ground installation. 


c. Cran iz IMO Or Tir PAPE 

The paper is divided into seven sections. This section has 
presented the problem and the objective of the study. Section II 
describes the elements of the TACTICS program and discusses the 
forces and the generation of these forces in the program. Section III 
compares the forces ona helicopter with the forces in the TACTICS 
peocram, describes the helicopter to be simulated im this study, 
and discusses the parameters and the derivation of these parameters 


used in the simulation. 








Basic maneuvers of the helicopter and the programing used to 
describe these maneuvers are presented in Section IV, along with 
several examples of the programing and computer output for the 
maneuvers. Section V presents two examples of the simulation of 
the hypothetical helicopter ina tactical situation, Section VI discusses 
modifications to the basic program, and the conclusions and recom- 


mendations are presented in Section VII. 








Il, THE TACTICS PRGGR 


The purpose of this section is to provide the reader witha 
general understanding of the concepts and theory of operation of the 
TACTICS program. This background will enable the potential user 
to comprehend more fully the difficulties involved in the simulation 
of a helicopter, and the reasoning behind the methods chosen. De- 
tailed information concerning the program is available in the TACTICS 
manual, Ref. 1, but it should be noted that TACTICS was designed to 
allow the analyst to obtain information without an in-depth knowledge 
ef the mechanics of the program. The TACTICS program can be 
considered to be composed of two elements, a group of computational 


routines and a Policy subroutine. 


x DE COMPULATIGCNAL ReV@TINGs 

The computational routines consist of the MAIN program and 
approximately fifty subroutines. The program reads the input data 
and initializes the starting positions and characteristics of three 
Iedies, normally a fighter, a missile, anda target. As time 1s 
incremented in the program, these vehicles will ''move'' in a straight 
line in three-dimensional space according to the imposed initial 
conditions unless acted upon by some external forces. These forces, 
namely aerodynamic, gravitational, and thrust, are specified or 
generated by the program, and through the familiar relationship 
F = ma, an acceleration is determined. The routines then integrate 
the equations of motion to determine the velocity and position of the 


vehicles and provide an output of the details desired by the analyst. 








The manner in which the forces act on the bodies is a result of the 
specifications in the Policy subroutine. 

The computational routines are not easily modified to any great 
extent, but as noted earlier, this element of the TAC [@saproc gar. 


does not have to be understood in detail to obtain useful results. 


LER Poole ew DROUIN 

The Policy subroutine is the means by which the analyst specifies 
the problem and dictates the desired performance of the vehicles 
under consideration. There are many guidance and control law ae 
tines availeble among the computational routines, and the Policy sub- 
routine is used to state the criteria for changing from one routine to 
mother. For exatnple, if the requirement exists for a vehiclesto ily 
straight and level for ten seconds, then start a 4.0 g climb, the 
sequence may easily be programmed in the Policy subroutine using a 
conditional GO TO statement based on time. Range and bearing 
angles are additional examples of criteria commonly used in the 
Policy subroutine. The applicable computational routines are listed 
in Appendix A, which includes a brief description of their function. 
When calling these subroutines in the Policy subroutine, the method 
in which the applied forces are to be developed is also specified. 
The development of these forces will be discussed in the following 


paragraphs. 


C, DETER Mi Tiel OT eenrers 
I. Phruse (1) 
There are three basic methods for determining the thrust 


force T ona body during any particular maneuver, the method 
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determined by the value of the variable dui. 

a. ITHR = 1 or 2 causes a variable value of military or after- 
burner thrust to be applied to the body, the value determined as a 
function of mach number and altitude through the use of tables. 

b. ITHR = 3 or 4 causes a constant military or afterburner 
thrust to be applied to the body, the value previously set in the 
initial data. 

c. ITHR = 5 causes no thrust to be applied. This mode is 
normally used when simulating a point mass. 

The value of thrust obtained by the above methods may be modified 
by the use of the variable THROTL. This parameter acts as a throttle 
setting and will adjust the thrust accordingly. The thrust force is 
assumed to act through the center of gravity of the body and along the 
longitudinal body axis. 

2. Drag (D) 

There are also three basic methods for determining the drag 
force D acting ona body, the method determined by the variable 
Eee RO. 

a. IAERO = 1 causes drag to be determined using analytic 


equations. The basic drag equations are: 


D eee 

Cine © eee Oceicient 

A = Aerodynamic reference area of the body 
- £4 

q = Dynamic pressure - 1b/ ft” 

9 = Air density - slug/ft” 

V = Magnitude of velocity of the body - ft/sec 


Ch is determined by the following expression: 


ee 











Ch = oe cf B Cr 
dic, ) 
C = Zero-lift drag coefficient 
D 

© 
ce A Re : 

D = A coefficient used with a parabolic 
acc, *) function for drag coefficient 


Menese coetiicients are read in as initial data. 


b. JAERO = 2 causes drag to be determined using the basic 
drag equations given above, but the value of the drag coefficient Cry 
poecetermined from tables as a function of mach number and lift 
fSecilicient. 

c. JAERO = 3 causes no drag on the body. This mode is 
normally used in the simulation of a constant speed vehicle. 


The drag force is assumed to act to directly oppose the velocity 


Srerie vehicle. 
Stee I) 
The lift force acting on a body is also determined using 


analytic functions. The following expressions are used: 








L = Cr DOUG 
dC 
Ce ese cence 
dg 
where Cr ite oOciiied ent 
dc, 
ee Slope of the lift vs. angle of attack curve 
a = Angle of attack 
% = Zero-lift angle of attack 


Zz 








The lift force is assumed to act through the center of etait. 
normal to the plane containing the velocity vector and the pitch axis 
of the vehicle. The angle of attack ( q ) is defined to be the angle 


between the thrust vector andytme velocityavector, 


4. Gravitational force 
In this problem, gravity is assumed to act normal to a flat 


earth, causing an acceleration of 32.174 cee 


Ns: 








Ii... DHE HELICOP Ta 


The TACTICS program was developed primarily for aerodynamic 
bodies having high speed capabilities. This section discusses the 
differences between the applied forces for these high-speed vehicles 
and the forces acting ona helicopter. A hypothetical helicopter is 
introduced to be used as an example, and the parameters used to 


describe it in the TACTICS program are presented. 


ee DISCUSSION OY PORCHES 

Since the TACTICS program was designed for the simulation of 
vehicles such as a fighter aircraft or a missile, the forces used in 
the program are specified and generated in the conventional manner 
explained in the following paragraphs. 

The thrust force acts along the positive longitudinal axis of the 
body to provide the desired acceleration along this axis. Conven- 
tional aircraft and missiles do not have the capability of changing 
Poe direction of the thrust in relation to the body. Deceleration of 
the vehicle along the longitudinal axis must be caused by a drag 
force greater than the thrust force, since the magnitude of this 
thrust cannot be made negative. 

The lift force acting on a conventional body is generated by an 
aerodynamic surface moving through an air mass. This force is 
considered to act normal to the velocity vector of the body and causes 
a lateral or turning acceleration. An aerodynamic surface, suchas 
a wing, having no relative velocity with respect to an air mass will 


Hou produce lit. 
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ee is a force produced by the body as it moves through the air 
mass and generates lift. It is considered to act basically along the 
longitudinal axis to oppose the forward motion of the body. As in the 
case of lift, drag is not produced by a body at rest in the air mass. 

The case of the forces acting on a helicopter during the hovering 
evolution, thatis, with the fuselage at rest with respect to the air 
mass, is now considered. Lift is produced by the action of the rotor 
blades moving through the air, and by proper adjustment of blade 
speed and angle of attack, the lift force is made to balance the force 
due to gravity, resulting ina net acceleration of zero. Within the 
confines of the TACTICS program however, lift is a function of the 
square of the velocity. This implies that as the velocity of the body 
approaches zero, the coefficient of lift must become infinite to sustain 
the Road in flight, an impossible condition. 

Vertical climbs and descents are an extension of the hovering 
evolution. The lift force must be varied to produce a net vertical 
acceleration, but with a forward velocity of zero, this is not possible. 

Investigation of force resolution in simple level flight produces 
similar problems. The helicopter may be considered to act asa 
body producing a single force vector. The component of the force 
vector along the flight path is labeled thrust, while the vertical com- 
ponent is considered to be the lift. Thrust and lift are therefore 
molated in a complicated way through blade speed, angle of attack, 
orientation of the rotor disc, and several other variables. Conversely, 
thrust is basically independent of lift in the TACTICS program. 
Problems of this nature also arise in relating drag forces on a heli- 


eempter to the drag force in the TACRICS prorvram. 


eS 








Consideration was given to the translation of helicopter forces 
to the forces in the simulation program. Attempts to isolate the 
generation of lift by using blade area and rotor speed proved unsuc- 
cessful however, as did direct control of the direction of the force 
vectors. Reference 3 discusses the simulation of.a helicopter rotor, 
but the method used was not considered practical for TACTICS 
applications. 

In nommal operation of the program, if the thrust force were 
set at a given value anda climb control-law routine were called, the 
TACTICS program would cause the vehicle to climb at the proper 
rate for the given thrust. Since it was not possible to resolve the 
helicopter forces and to determine the proper coefficients, the ap- 
proach taken was to simulate the helicopter by specifying the perfor- 
mance parameters, rather than the forces. For example, this 
method specifies the proper rate of climb for the helicopter under 
consideration using published performance curves, then the appropri- 
ae oOrces to produce this rategot climb are developed. in effect, 2 
fixed-wing vehicle is used having theoretical parameters that enable 
Metoeperform like a helicopter. 

In general, the motion of the body is controlled by proper ad- 
justment of the thrust force to obtain a desired acceleration, and by 
Sontrol of the velocity vector for the correct direction. Specific 


maneuvers and their operation are discussed in detail in Section V. 


I, HELICOPTER CHARACTERISTICS 
Reference 4 discusses methods for the performance analysis of 


an existing or contemplated helicopter. Reference 5 provides a 
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sample performance calculation for a theoretical helicopter using the 
guidelines of Ref. 4. While the performance capabilities of this 
helicopter are somewhat less than those of a modern helicopter, it 

is felt that the methods of analysis used are applicable to the study of 
the newer versions. For this reason, the hypothetical helicopter of 
Ref. 5 was chosen as the demonstration model for this simulation. 


The pertinent characteristics of this vehicle are shown below. 


opecifications for the Hypothetical Relveopren 


Type Single-Rotor Cargo Utility 


Rotor Diameter 52.0 feet 
Gross Weight Range 5000-7000 pounds 


A summary of the performance characteristics for a gross 


weight of 6300 pounds is given below: 


Vertical rate of climb at sea level (20 Thome 
Hovering ceiling (Out of ground effect) 5500 feet 
Best climbing speed at sea level 55 mph 
Best rate of climb at sea level 1450 fpm 
Service ceiling 13,300 feet 
Maximum level flight speed at sea level 118 mph 
Average cruising speed 84 mph 


Since the basic system of measurement in the TACTICS program is 
the foot-pound-second, some of the specified values must be con- 


verted to this system. 
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ce: PARAMETER DEFINITION 
This section contains a description of the parameters used in 
the TACTICS program that cause vehicle three to assume the 
characteristics of the hypothetical helicopter. 
1. DATA 44: W0(3) = 6000.0 
Initial weight was set at the average value of the gross weight 
range. Since fuel burning rate was not considered in the simulation, 
this value remained constant. 
CeeeD ATA 49 meR EA Giese 
The aerodynamic reference area was arbitrarily selected to 
be that of the rotor disc. Since the program assumes a given perfor- 
mance capability, rather than calculates it, this value is not critical. 
If the reference area is selected to be a very small value however, 
the lift coefficient must correspondingly become large to sustain the 
vehicle in flight. This can cause overflow to occur during execution 
of the program since the value of the lift coefficient is squared in 
several program statements. 
Ore oe 4 ee vi es)) =e 0 
Maximum lateral acceleration of the theoretical helicopter 
was estimated to be 3.0 g's. 
ae Uae Oe Vie ) a 10° 
The maximum coefficient of lift was assumed to be infinite 
eg. ; 10°). This was necessary in order for vehicle three to develop 
pene required lift at low airspeeds. 
5. DATA 91: CDOCON(3) = .014 
It is difficult to directly compare the drag forces produced by 


eeneINcOpte rama cd lime wine alle rattinatiec sa mic irelere nee ml raiiic. 
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A fixed wing aircraft at rest in the air mass obviously does not produce 
drag, nor does it fly. The rotor blades of a helicopter ina hover, 
however, produce profile drag and induced drag while generating lift. 
At higher speeds, both bodies are subjected to parasite drag, caused 
by the friction of the air mass on the fuselage. This element of total 
drag is the predominant force opposing motion of the vehicle in flight, 
and it is possible to directly translate the parasite drag of the hypo- 
anetmcal helicopter into tae LAC TICS meterencestrame: 


Table VI of Ref. 5 is produced in part below: 


v RHP 
59 12 
88 39 

117 88 

147 168 

176 286 

205 444 


where V is the velocity of the helicopter in feet per second, and RHP 
is the horsepower required to overcome parasite drag. To convert 


the power to actual drag in pounds, we use the equation 


SOOPER 


S V 


Using the analytic functions, drag is determined by the expression 


Dek C.kp Wo 


Simultaneous solution of these equations yields 


ee comes 


D ey 
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Using the values listed above for RHP and V, A equal to AREFA(3), 
and an average value of 9 equal to . 002, Ch is found to range only 
from .013 to .015. An average value of .014 was selected for this 
dima’ cocifeicenn: 
6. DATA 92: BCON(3) = 107° 

This 1s a coetticient relating the change an thegd ia cscoer 
ficient to the square of the lift coefficient. It may also be used to 
relate the increase in drag caused by the increase in angle of attack 
as a fixed wing body decreases in velocity. Since the induced drag 
tmot considered for the helicopter, this cociticienticmomd be zero, 
During program execution, however, this coefficient appears in the 
denominator of certain expressions, and was therefore set ata 
small number other than zero. 

7. DATA 93: SLOPE(3) = 20,000 

DPhistis a -coctiicient relating the Change iii eocmlcicn: 
with a change in angle of attack. Fora given airfoil, this is a 
meotimite value, oimece the liit coefficient for the helicopter hasva 
range that theoretically becomes infinite, the angle of attack must 
also become infinite to achieve this lift coefficient. Therefore, the 


large value of ae was selected to prevent the angle of attack 


da 
from becoming large, and has the effect of maintaining a level fuse- 
lage attitude even at very low velocities. 
8. DATA 94: ALPHAO(3) = 0.0 


Zero-lift angle of attack was arbitrarily set at zero, that 


Tope vel gicmt: 


20 








9. DATAI05: “THCONG)==. 800070 
Since thrust must exceed the weichtteia ped, an ond cimien. 
the body to accelerate in a vertical climb, eonstamt mm lita tae 
was set at a value to provide 2000 lbs. excess thrust if desired. 
Since the acceleration is a limiting value in this simulation, and the 
thrust required to produce this acceleration is obtained through the 


ase Of the throttle parameter, this value is not critical, 


a 








IVS HE LICOPLER MANB Nias 


This section describes the programing used to accomplish the 
simulation of the basic maneuvers of the hypothetical helicopter. It 
was the basic objective during the writing of these maneuver routines 
to test the feasibility of the method under consideration. Since it was 
extremely difficult to trace a given sequence of operations through 
the computational routines, it was generally not known in advance if 
a particular segment of programing would produce a desired result. 
During this debugging phase, the computational routines were com- 
piled and placed in permanent storage in the computer system. All 
programing to test a given maneuver was done in the Policy sub- 
routine, which was then compiled separately, linked to the computa- 
tional routines, and executed. The use of the Policy subroutine in 
this manner was justified by the ease with which a programing seg- 
ment could be modified, and by the lack of detail actually desired. 
An expanded and more detailed version of the helicopter program 
would probably use individual subroutines to accomplish a given 
maneuver, Similar to the control-law subroutines incorporated for 


Peebter aircratit. 


ee PAE HOVER 

To preclude any difficulties in the generation of lift at zero 
airspeed, the helicopter is required to maintain a velocity greater 
than the arbitrarily set value ef 0.5 ft/sec. It is felt that for all 
practical purposes, this velocity has no effect on a tactical situation, 


and the helicopter is essentially ina hover. Once this speed is 


oe 








attained, it is possible to maintain the hover by balancing the 
gravitational force with an equal and opposite thrust force. Since 

the vehicle is essentially a motionless inass in three-dimensional 
space, however, itis treatedas such. By calling the subroutine 
STRFLT(3, AERO, ITHR), with AERO = 3 and ITHR = 5, the vehicle 
will maintain a constant direction with no change in velocity. 

It is of interest here to discuss the effects of certain sub- 
routines when the vehicle is ina hovering condition described above 
or at a low velocity. Suppose the vehicle is in a hover with the 
velocity vector ina horizontal plane with a magnitude of 0.5 ft/sec. 
It is desired to have the vehicle climb at a 45 degree angle from the 
hover, requiring that the velocity vector be rotated 45 degrees above 
the horizontal. This can easily be accomplished by using the sub- 
noutine CLIMBI(I, GFORC, IAERO, ITHR), which is a constant- 
velocity, constant-g climb, but consideration must be given to the 
argument GFORC, the amount of g specified. A value of GFORC 
equal £0 3.0 will almost instamtly cause the velocityiyvector to exceed 
the 45 degree angle desired, since there is a large force (three 
times the body weight) applied to the body in the vertical direction. 
To prevent this large overshoot, GFORC should be set at a value 
near to 1.0. If diving and turning subroutines are used at low 
melocities, Similar problems occur. Using values near unity permits 
the use of conditional IF statements in the Policy subroutine to 
monitor the velocity vector angle and prevents large discrepancies 
between desired and actual angles. The following discussion of 
vertical flight has an example that terminates with the hypothetical 


Helveeprer imine vie 2: 
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i, VERTICAL FLIGHT 

To perform vertical maneuvers, the vehicle imjtehe, wae mie. 
program must be oriented with the thrust vector parallel or anti- 
parallel to the gravitational force vector. Using analytic functions, 
some lift will be generated as the body ascends or descends, causing 
aenorizontal component of the resultant foree vector scr mp iemone 
body. This component will be small due to the low velocities associa- 
ted with vertical flight, and can be neglected or considered a minor 
error inthe pilot's ability to maintain perfectly vertical flight. 
pwmeuld the meed arise for completely vertical flight, the velocity 
vector can be adjusted to provide a component of thrust to oppose the 
lift, but in this initial study, this accuracy was not justified. 

Reference 5 calculates the vertical rate of climb at sea level 
for different gross weights. For the theoretical helicopter, this 
value is 1150 ft/min at 6000 pounds. Although calculations for the 
vertical rate of climb as a function of altitude were not performed, 
Ref. 5 indicates the general shape of the graph, and Figure 14-7 of 
Ref. 2 shows the curve for anactual helicopter. It may be noted 
that for gross weights in the mid range and altitudes below 7000-8000 
feet, the rate of climb is essentially constant. 

In the simulation, the variable VERCLB was set to 1150.0 and 
Hope transtorimed to fect per second, Since the vertical climb is 
mee d mainly for takeoffs, the constant assumption in this simulation 
ieecousidered appropriate. = Por more advanced use, it is a relatively 
simple matter to make a table of VERCLB vs. altitude and thus 
determine the rate of climb appropriate for any given altitude. This 


technique is demonstrated in the section covering autorotation and 


forward climb. 
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Since the ability to increase the velocity im @ yerticaledirectionm 
is dependent upon the thrust over and above that thrust required to 
counteract the body weight, it was decided to use vertical acceleration 
(caused by this excess thrust) as a parameter in determining the 
vehicle's vertical movement. Although this data was not available 
for the hypothetical helicopter, it is assumed that this data can be 
calculated or determined through flight tests, and an estimated value 
oe, 5 cee was used to initialize the variable VACCEL. Similarly, 
a variable VDECEL was initialized at 27.0 Rieece LOmindieatesune 
mia cima rate of acceleration ior a verticamy descending body.) It 
was assuined that the helicopter would not be permitted to reach the 
ec ic tronic at cae toiecntroleree Er ienone) 

While the variables described above were initialized in the 
Eolecy subroutine, it would be a relatively simple procedure ¢o 
modify the computational routines to allow the accelerations and 
climb rates to be read in the same form as the standard input. 
Acceleration capabilities could also be related to altitude if this 
mecuracy 1s desired. 

Once these limits on vehicle performance have been established, 
it is necessary to determine the thrust required to execute the ver- 
tical maneuvers. The force required to accelerate a mass M with 
meee leravion VACCE lis M “= VACCEL. Nor the helicopter, 

M = WEIGHT(3)/G, and therefore the net force required for vertical 
acceleration is WEIGHT(3)*VACCEL/G. Neglecting lift, this net 
force is the resultant of the thrust, drag, and gravitational forces 
acting along the vertical flight path. Modification of equation (10) 


of Ref. 1 yields 


Zo 








WEIGHT(3)*VACCEL/G = T * COS a - DRAG(3) - WEIGHT(3)*SIN Y. 


Since IT = THROLEL(3)*1 HCON(3), the throttle setting neecessan are 


cause an acceleration of VACCEL is determined t,, 


WEIGHT (3)*VACCEL 


eNom (3h) = G + DRAG(B)AW FIGHT SUNG (ao 


—_— —= 


THGOM(s) COs (Ar eri wy) 
Similarly, the throttle setting for a constant rate of climb is de- 
termined by setting thrust equal to the drag and weight as follows: 


- WEIGHT (3)*SIN(V(3, 6)) + DRAG(3) 


inmonder onde celétatre atsaeuatrce.y DiKereeathe throttle pa yamieter 


must be decreased to the setting 


= WaetIGHT(3)*VDE CEL 
THROTL(3) = G + DRAG(3)4W BIGHT (3)-SINiu. ca 


———— 


THCON(3)*COS(A LPHA(3)) 


The statements listed in Fig. 1] pertain to the hypothetical heli- 
copter and have been extracted froma Policy subroutine. These 
statements will cause the helicopter to: 


a) Perform a vertical takeoff after one second of elapsed 
time 


b) Accelerate at VACCEL eee until maximum vertical 
climb rate of VERCLB ft/min is attained 


c) Climb vertically to an altitude of 85.0 feet 


d) Decelerate at VDECEL fee to a hover 


Initial conditions established the vehicle at an altitude of one foot with 
a horizontal velocity of 0.5 ft/sec. Fighter and missile values were 
Ectiat ZeLOs ss CO puLe mG plOc Ol tae =nosition oe! the Meliconve mn imitie 


vertical plane is shown in Fig. 2. The position was plotted at 0.5 
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Extract from Policy subroutine with programing 
for vertical climb maneuvers. 
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Figure l., 
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second intervals, with the unequal spacing of the points indicating 


accelerations. 


c- HORIZON iA bh LIiGH & 

Horizontal flight is accomplished by a method analogous to that 
of vertical flight. Inthis case, a horizontal acceleration and de- 
celeration capability of HACCEL = 8.0 eee and HDECEL = 5.0 
aera is assumed. The vehicle is oriented with the velocity vector 
in the horizontal plane and the throttle is adjusted to produce the 
thrust required for the desired acceleration. In this case, drap be - 
comes more noticeable due to the higher velocities attained, and 
acceleration will be somewhat less than in the absence of drag, white 
deceleration will be correspondingly higher. The thrust force in this 
case is assumed to be a constant, dependent on the position of the 
rotor disc. Since the rotor disc may be tilted backward, the net 
force produced by the disc may have a component directly opposing 
forward flight. This condition does not exist wititfixed wing aircrait, 
but it can be simulated in the program by the use of a negative throttle 
parameter. 

A helicopter in forward flight has many of the characteristics 
of a conventional aircraft, and can be simulated in the TACTICS 
program using the normal analytic functions. For example, a diving 
turn can be accomplished by lowering the velocity vector below the 
horizontal using a DIVE subroutine, then calling a TURN subroutine 
to complete the maneuver. As before, acceleration along the flight 
path must be controlled by the thrust, unless constant-speed routines 


are used. Excess thrust during a diving turn will cause the vehicle 
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to accelerate, and if the turn is specified as constant-c, thesmrahn 
path will describe amincreasing spiral. 

Use of a TURN subroutine when the vehicle is at a very low 
velocity has the effect of pivoting the vehicle about the vertical axis. 
A helicopter can easily accomplish this maneuver, and it may prove 
useful in providing a stable platform from which an opponent can be 
continuously tracked. 

The extract from a Policy subroutine listed in Fig. 3 shows the 
statements that will cause the helicopter to: 

a) Decelerate to a hover at HDECEL foe 
b) Hover for ten seconds 
c) Accelerate at HACCEL fy ceca 


d) Start a diving right turn of 1.2 g's after attaining a 
velocity of 88.0 ft/sec (60 mph) 


Initial conditions located the helicopter at 4000 feet with a velocity of 
60.0 ft/sec. A computer plot of this maneuver in the horizontal plane 
is shown in Fig. 4. The apparent left turn of the vehicle is due to the 
orientation of the axes in the TACTICS program, as explained in 

Ref. 1. As before, the uneven spacing of the position points plotted 


at one second intervals indicates accelerations of the vehicle. 


D. FORWARD FLIGHT CLIMB 

For a given gross weight and optimum climb speed, the rate of 
climb of a helicopter will decrease with altitude. The service ceiling 
of the helicopter is determined by that altitude where the rate of 
climb isless than 100 ft/min. This section demonstrates a method 
that includes a table of rate of climb vs. altitude to simulate a realis- 


tic forward climb flight path. This method may be applied to other 
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maneuvers in the simulation, suchas the vertical rate of climb vs. 
altitude problem mentioned earlier. 

Figure 19 of Ref. 5 illustrates forward rate-of-climb capability 
vs. altitude at normal rated power for the hypothetical helicopter. 
This graph was used to initilize an array named FWDROC with paired 
values of rate of climb and altitude. During program execution, the 
array is scanned to find the proper rate of climb for the current 
ercitude of the helicopter, The velocity vector of the vehicle is then 
adjusted by CLIMB or DIVE subroutines to make the vertical com- 
moment of this vector equal to the proper rate of climb. Since there 
is normally some overshoot due to the integration step size when 
adjusting to the desired climb angle, an error of up to 24. 0 ft/min is 
allowed in order to prevent constant climbing and diving about the 
eorrect value. 

Optimum forward velocity for the best rate of climb is deter- 
mined by constructing a horizontal tangent to the power-required 
curve for the hypothetical helicopter, Fig. 10 of Ref. 5. The best 
climb speed determined by this technique is 55 mph. The helicopter 
is normally accelerated to this velocity in level flight prior to 
faitiating the climb. 

An example of this climb method is illustrated in Fig. 5. The 
iaitial position of the helicopter is at 5000 feet of altitude with a 
locity of 40.0 ft/sec. Figure 6 is a plot of the Aight path of the 
vehicle in the vertical plane, and the decreasing slope of the curve 
indicating a decreasing rate of climb may be noted. The higher 
altitudes were chosen for the example, since the climb rate is almost 


constant in the lower ranges. 
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This method is considered adequate for the purposes of this 
paper. For amore detailed analysis, interpolation routines should 
be used to provide a functional relationship between rate of climb 
and altitude, and the process would more appropriately be located in 


a subroutine. 


i. BOR VAR bien AU TORO TATICGnN 

Autorotation of the hypothetical helicopter is simulated by a 
method analogous to that employed in the forward climb maneuver, 
Put in this case, the rate of descent is a finetion of forward aoe 
Figure 15 of Ref. 5 was used to initialize an array named AUTORO 
with paired values of rate of descent and forward speed. This array 
was then Scanned to determine the rate of descent corresponding fo 
the actual horizontal component of velocity of the helicopter. The 
angle of descent is then adjusted by CLIMB and DIVE subroutines, 
and the thrust regulated to cause the vehicle to follow the correct 
flight path for an autorotating helicopter. A more realistic approach 
to the programing of the method during an analysis would be to use 
subroutines and interpolation routines for a more accurate approxi- 
mation of the flight path. 

An example of the programing used to simulate the autorotation 
maneuver is used in the Policy subroutine demonstrating an auto- 
rotating helicopter under attack by a ground-to-air missile, Appendix 
Bb, 2ection Hl A computer plot of the maneuver is given in Section 
IV of the same Appendix, but the scaling factor causes the maneuver 


to appear steeper than it actually is. 
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V> TACTICAL PE XAMPIEBS 


Two examples indicating potential applications for the methods 
developed in the simulation of a helicopter are presented in this 
section. The examples illustrate Policy subroutines that use com- 
binations of the maneuvers described earlier, but no attempt was 
made to determine optimum maneuvering for either vehicle. An 
analysis of a tactical situation would normally consider only certain 
aspects of a given attack. Later developments would be considered 
separately. The fighter vs. helicopter example is extended beyond 


the initial attack by the fighter for demonstration purposes only. 


i . DiGi he Vo. EIR TE 

In this example, a fighter aircraft with an altitude advantage 
makes a diving gunnery run on the hypothetical helicopter, then pulls 
up and circles around for another attack with 20-mm projectiles. 
The fighter aircraft simulated is an F - 104, and tables are used to 
describe the interactions among the lift coefficient, drag coefficient, 
angle of A Ae. maximum coefficient of lift, and mach number. The 
Policy subroutine defines the problem as follows: 

VERICLE 1(F - 104¢pumsuine aircraft) 


- Fly pursuit course navigation against helicopter; 
throttle set at .9. 


- When relative range to helicopter is less than 4500 ft., 
fire two 20-mm projectiles at the target. 


= litvaltitude falls belew 3500 it... “climb at Sa5e2 until 


nose is 3 degrees above the horizon; throttle at full 
military power. 
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- Performa 4.5 g climbing right turn until bearing to the 
helicopter is less than 60 degrees, then fly pursuit course 
navigation with throttle set at . 80. 

- When relative range to helicopter is less than 4500 ft., 
fire two 20-mm projectiles and climb out at 5.5 g when 
altitude falls below 1000 ft. 

SC (20mm projectile ) 

- Fly captive flight until launch criteria are satisfied. 


- Launch and fly in a ballistic trajectory. 


- Find miss distance, restore conditions to launch time, 
launch second projectile 0.01 seconds later. 


VEHICLE 3 (Hypothetical helicopter) 


- Perform constant 1. 3-g turn toward the fighter to 
attempt to force a head-on passing situation. 


- When fighter passes and starts to turn for re-attack, 
slow to 25 mph if possible. 


- When fighter begins second tracking run, autorotate at 
current forward speed. 


The initial position of the fighter is at the origin at an altitude 
of 7000 ft., with a heading of 090 degrees at 800 ft/sec. The heli- 
eopter is located 2,000 it. directly East at an altitude of 3000uC)| 
heading 350 degrees at 120 mph. The complete Policy subroutine 
Zor anes tracted moolein 1s liste deine sppenaiy.-) oe Obl Om) me Ehime 
orthogonal views of the flight paths of the fighter and the helicopter 
are shown in Section II of Appendix B. The points were plotted at 
evo second auntervals by the computer, and the connecting line 
segments were then added to make visualization of the flight paths 
easier. The turning helicopter causes the fighter to make constant 
corrections during the initial attack, and since the fighter is using 
pure pursuit techniques, rather than adding some type of lead cal- 


culation, the miss distance of the 20-mm projectile is approximately 
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100 feet. Similar results are obtained during {the secendmanac. 
when the helicopter is autorotating at approximately 40 mph. The 
attacking figher was ina position for re-attack before the helicopter 


could slow all the way to 25 mph. 


B. MiSSIpE ey >. tie LICOrrER 

This example demonstrates the capability of the TACTICS 
program to simulate the helicopter in a ground-to-air missile envir- 
onment. Due to the high velocity of the missile and its proportional 
navigation system, the helicopter attempts to evade the threat by 
autorotating at high speed when the relative range is less than 2.5 
miles. The helicopter is initially 25, 000 ft. away from the missile 
site at 1000 ft. of altitude. It is heading directly at the site ata 
speed of 120 mph, and at missile firing, speed is increased as much 
as possible within a 140 mph maximum prior to commencing the 
autorotation. The ground-to-air missile has the characteristics of 
the vehicle described in Example 3, Ref. 1. The Policy subroutine 


defines the problem as follows: 


VEHICLE 1 (Not used) 
V DPIC LE 2. ou rtace= Osan missile} 


Pp igawncm at times zero; thrust constant at 5577 lbs] initian 
velocity of 100 ft/sec. 


- Fly in O-g dive until burnout at 1.0 seconds. 
- Continue O-g until time to guide; no thrust. 


- Begin guidance 1.0 sec. after launch; fly proportional 
navigation with time lag of 0.2 sec; no thrust. 


VEHICLE 3 (Hypothetical helicopter) 


- Fly straight and level at missile site until missile launch. 


oo 








- Accelerate at maximum capability until range to missile 
LS) C.g5 eines. 


- Autorotate at current forward speed. 


The Policy subroutine for this scenario is given in Appendix B, 
Section III. The computer plot of the flight paths of the two vehicles 
in the vertical plane is shown in Section IV of Appendix B. The heli- 
copter was able to increase its velocity only four mph prior to 
initiating autorotation. As expected, the miss distance of the missile 
was only 1.5 feet, due to the slow speed of the helicopter and the 


proportional navigation system used by the missile. 
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Vi. PROGRAM MODIFICATIONS 


Minor modifications to the TACTICS MAIN program were made 
to facilitate input/output operations on the computer. Reference 6 
discusses operation of the TACTICS program at the Naval Postgraduate 
School, and was used as a reference for program execution during this 
study. A subroutine used at the Postgraduate School was also modified 
and added to the computational routines to add a plotting capability to 


thiesoroeram, these modifications will be discussed im this Section: 


as MAIN PROGRAM 

The IBM 360/67 computer system at the Naval Postgraduate 
School was utilized for the execution of the program. The system 
normally operated under two modes, the TSS/360 time sharing system 
and the basic OS/360 system. Since the input and output statements 
differed in each mode, the variables INGAGE and OTGAGE were used 
in the READ and WRITE statements of the MAIN program to facilitate 
the change from one mode of operation to the other. These variables 
were also used in each subroutine having READ and WRITE state- 
ments, and were placed in COMMON/GAGE/. Under TSS/360 
operation, the MAIN program asks for the values of INGAGE and 
@eGAGE, and the operator enters these values at the remote ter- 
minal. Under OS/360 Gpe tation meierdcoared va liioscm ese ite hearas 
ene first data card. 

The main program was also modified by the addition of state- 
ments calling the subroutine OSPLOT at the termination of vehicle 


maneuvering. The arrays used as calling arguments in OSPLOT 
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were placed in COMMON/PLOTT/ and initialized to the value 900, 000. 0 
at the beginning of the MAIN program. An explanation of these arrays 
and the function of OSPLOT is given in the following paragraphs. The 


modified MAIN program is listed in Appendix Cayce cio. 


B, SUBROU TiN, OSEeOT 

OSPILOT is a subroutine that was written at the Naval Post- 
graduate School for the purpose of plotting sets of points ona single 
graph. It was decided to use this subroutine to provide the user with 
a rough plot of the flight path of the vehicles under Pon eeneiom 
and to permit visualization of the maneuvers without scanning pages 
of computer output. Since the variable DIT PO determines the time 
interval between information outputs, it was a logical choice for the 
time spacing of points on the flight path plot. During program execu- 
tion, the Policy subroutine stores the coordinates of the vehicles, 
auethe chosen points in time, in the arrays of COMMON/PLOTT/. 
This information is then used by OSPLOT at the completion of man- 
euvering to produce a two-dimensional plot of the flight path. The 
subroutine is normally called three times in the MAIN program with 
the arguments arranged to produce three orthogonal views of the 
maneuvers. 

During the experimentation phase of maneuver testing, the 
automatic scaling features of OSPLOT were used for ease of operation. 
When further refinement of the plot was desired, manual scaling or 
adjustment of DT PO was employed. The printing of the axes was 
eliminated from the subroutine to make the output neater. 

The subroutine was also modified to print different characters 


for each vehicle. OSPLOT requires that if the abscissa and the 
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ordinate of one vehicle are known and are to be plotted, the ordinate 
corresponding to the same abscissa for the other vehicle must also 
be known. Since the vehicles normally have no points in common, the 
value 900, 000. 0 (initialized in the MAIN program) was given as the 
ordinate of the second vehicle. When the program scans the co- 
ordinates for plotting, values greater than 800, 000. 0 are ignored, 
ene the ficthecrous points are not plotted. The revisediversion of 


@=ePLOT is listed in Appendix C, Section II. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


The primary difficulty in applying the TACTICS prepraine tome 
simulation of a helicopter in a tactical environment is the translation 
of the forces produced by the rotating blades of the helicopter to the 
conventional fixed-wing forces used in the TACTICS simulation pro- 
gram. After a study of this problem, it was concluded that ifa 
transformation of the forces was even possible, it would require 
major new subroutines to be added to the existing Ue areal. 
routines. It was therefore decided to use the known capabilities of 
the helicopter as limiting values, and te use the computational rou- 
tines already available to cause the vehicle to perform within these 
limits. This method utilizes the capabilities of the TACTICS program 
and allows helicopter and fixed-wing vehicles to be ''flown" in the 
Same simulation. 

Additional work remains to be done in the programing of heli- 
copter maneuvers, and the flexibility of the program in the simulation 
wimtie, Welicopter can be improved, The use of excess thrust in place 
of acceleration as a limiting factor is another area requiring additional 
mivyestigation. Nevertheless, the ideas presented in this paper are 
feasible and represent one approach to the problem of simulating a 


helicopter in a tactical environment using the TACTICS program. 
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Pol eis ND Ex 7 


TACTICS SUBROUtIVE DESCRIPTIONS 


Subroutine Desicripion 
MAIN Zeros values in COMMON/TABLE/ 


Determines number of cases in run 
Determines location of input and 
output files; times run 


INCOND Prints initial data values 
Sets initial flays and zeros out 
variables 
Computes relative range rates and 
VElLOCity VECciors 


So CPU TL Prints output as controlled by NPRINT 
options specified in POLICY 


AUXCOM Computes auxiliary computations for 
printing only 


ieee lb R Reads tables when provided for in 
vehicle's aerodynamic characteristic 
description 


DE CRD Reads data cards for information as 
provided on input data form 


Pe IT UD Computes attitude angles 

CLIMB1 Elementary climb where G-force is 
specified 

CLIMB2 Constant speed climb ; depends on 
excess thrust 

DIVE] Simple dive, G-force is specified 

SORE LT Vehicle flies in a straight line but 
not necessarily level 

STRLVL Vehicle flies straight and at a constant 
alitiudcde 

RTRNI Simple right turn; G-force specified 


determines angle of bank 


RTRNZ Constant-altitude, constant-Mach right 
turn; depends on excess thrust 
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Subroutine 


Rebs 
LTRNI 
LTRN2 
LTRN3 
AIM 
TABINT 
INTGRT 
LIMIT 


eA PFLT 


ATMOS 
WEIGH 
THRUST 


LAG 


AERODN 


INITS 


PLACRD 


Des CEiption 


Constant- Mach right turn where G-force 
is specified 


Simple left turn; G-force specified 
determines angle of bank 


Constant-altitude, constant —-Wwiaci lee 
turn; depends on excess thrust 


Constant- Mach left turn where G-force 
is specified 


Aims vehicle one as specified in initial 
conditions 


Interpolates values as provided in 
input data tabies 


Computes relative and absolute positions, 
velocities, and ranges 


Computes aerodynamic limitations as 
specified in the input data 


Zeros out variables for vehicle(s) not 
used or when attached to another vehicle, 
such as a missile 


Computes model atmosphere, speed of 
sound and Mach number for given altitude 


Computes weight of vehicle as fuel is 
consumed 


Computes thrust available either from data 
tables or an input constant 


Computes change from commanded lateral 
acceleration (GFORC) to actual output 
acceleration 


Computes maneuver-induced drag, input 
acceleration components, and changes 
in velocity due to maneuvers 


Initializes conditions for ARKAM 
inte gration 


Computes throttle setting based on air- 


Grait limitations, | olicy 1npuL, and 
atmospheric conditions 
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Subroutine 


SHO 


FUNCTN 


MAG 


COORD 


DA UX 


ARKAM 


MA CHR8& 


rOPCON 


WOT 


RATES 


CROSS 


BZOMM 


in ONAV 


PR UIT 


LAUNCH 


esPLlor 


Desc ription 


Stores position and velocity when using 
restore option 


Computes unit vector components of com- 
manded acceleration and thrust 


Computes magnitude of lateral com- 
manded acceleration 


Makes transformations from rectangular 
to spherical coordinates and vice-versa 


Integrates cartesian and geocentric 
coordinates for all vehicles relative 
LOM point of intercept 


Integration subroutine 


Calculates aerodynamics in reverse for 
constant- Mach maneuvers 


Computes local and inertial coordinate 
transformations 


Computes vector dot products 


Computes range and angular rates of 
change 


Computes vector cross products 


Computes ballistic trajectory of 20-mm 
projectile 


Proportional navigation subroutine for 
intercept problem 


Pursuit course navigation subroutine 
for intercept problem 


Subroutine called for the launching of 
a missile or projectile from vehicle 


Subroutine that outputs two dimensional 
plot of vehicle positions 
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